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Major observational efforts in the coming decade are designed to probe the equation of state of dark 
energy. Measuring a deviation of the equation-of-state parameter w from -1 would indicate a dark 
energy that cannot be represented solely by a cosmological constant. While it is commonly assumed 
that any implied modification to the ACDM model amounts to the addition of new dynamical fields, 
we propose here a framework for investigating whether or not such new fields are required when 
cosmological observations are combined with a set of minimal assumptions about the nature of 
gravitational physics. In our approach, we treat the additional degrees of freedom as perturbatively 
constrained and calculate a number of observable quantities, such as the Hubble expansion rate 
and the cosmic acceleration, for a homogeneous Universe. We show that current observations place 
our Universe within the perturbative validity of our framework and allow for the presence of non- 
dynamical gravitational degrees of freedom at cosmological scales. 



I. INTRODUCTION 

Cosmological observations suggest the existence of 
dark matter. There is, also, by now ample evidence that 
the Universe at late times is accelerating. Do either of 
these facts require the introduction of new, dynamical 
degrees of freedom in the gravitational and matter sec- 
tors? 

In the case of dark matter, the analysis of "bullet clus- 
ter" objects such as IE 0657-56 Q and MACSJ0025.4- 
1222 comes as conclusive evidence for the existence of 
new degrees of freedom independent of the visible bary- 
onic matter. While the exact nature of these new degrees 
of freedom are still under debate, the simplest choice - a 
single, non-relativistic, dark matter fluid - is the center 
of the cosmological standard model. 

Even if we choose to account for the galactic rotation 
curves with a modification to gravity, it is difficult to 
do it without the addition of extra degrees of freedom. 
A consequence of the Lovelock-Grigore theorem [3|, |4J is 
that most local, classical, and covariant modifications to 
the gravitational action lead to new fields [|. Alter na- 
tives to dark matter that have these properties - such as 
TeVeS [1] and STVG Q - themselves have new fields as 
a consequence of the Lovelock-Grigore theorem. 

Today, cold dark matter with a cosmological con- 
stant (ACDM) satisfies all current observational con- 
straints [I]. However, a number of theoretical issues @ 
have prompted observational campaigns to search for de- 



viations from the predicted ACDM behavior [T(|. These 
observations will measure the equation-of-state parame- 
ter w of the dark energy and search for deviations from 
the value w = — 1 that corresponds to a cosmological con- 
stant. Should their results force us to abandon ACDM, 
we will face a problem for dark energy similar to that we 
faced for dark matter: do such observations require new 
dynamical degrees of freedom? So far, the two major ap- 
proaches to dark energy modeling that predict values of 
w 7^ —1, quintessence [III an d f{R) gravity pjj], incor- 
porate new dynamical degrees of freedom, which can be 
cast, in both cases, in the form of a single scalar field. 

A number of attempts have been made for modify- 
ing General Relativity in a way that generates cosmolo- 
gies with but without introducing new dynam- 
ical fields. Modified source gravity [13J, based on in- 
sights from the Palatini formulation of the Einstein equa- 
tions [14| forces the (Einstein-frame) scalar field to be 
non-dynamical by erasing the kinetic term. Cuscuton 
cosmology [TH, [l6[ also modifies the kinetic term; it ap- 
pears to have new dynamical degrees of freedom but these 
are frozen out at the perturbative level. Non-local p7l[T8j 
or holographic theories may also avoid the introduc- 
tion of new degrees of freedom. Albeit ambitious, these 
attempts can face serious problems (20| . 

In the opposite extreme, because of the Lovelock- 
Grigore theorem, the theoretical pressures on purely phe- 
nomenological models that deviate from ACDM but that 
do not introduce new dynamical degrees of freedom are 
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strong. The fact that observations favoring dark energy 
are made on relativistic, Hubble scales makes it even 
harder to avoid the consequence of the Lovelock-Grigore 
theorem, while remaining internally consistent. Requir- 
ing general covariance from the start, for example, pre- 
vents one from even phrasing many minimal theories the 
way MOND did on galactic scales [2l[ . In general, while 
many theories similar or equivalent to Brans-Dicke [33j 
may, in the limit of small post-Newtonian corrections, 
appear only to modify the General Relativistic equations 
of motion, the new scalars such theories introduce appear 
as truly independent degrees of freedom in the relativistic 
regime - for example, on scales approaching the horizon 
size. 

Approaches, however, that tie together two particu- 
lar sets of observations - homogeneous expansion and 
linear structure formation - in a phenomenological but 
self-consistent fashion have also been recently explored. 
For example, the parametrized post-Friedmann frame- 
work [22I makes minimal assumptions about the under- 
lying physics, requiring only causality, metric structure, 
and e nerg y- momentum conservation. The Cardassian 
model [23} introduces a power-law term in the Friedmann 
equation, which leads to accelerated expansion even in a 
matter dominated universe. The approach of Ref. [fj is 
similar in spirit, and explicitly shows how a subset of 
modifications will introduce new degrees of freedom. Al- 
though these approaches are general and simple to use in 
connection to observational data, they can not, by their 
very nature, generate predictions for other astrophysi- 
cal settings - the solar system, for example, or nearby 
(non-cosmological) event horizons. As a result, they can 
not guarantee that dynamical degrees of freedom are not 
required for consistency in these other situations. 

We follow here a different approach intermediate to the 
two outlined above, which we call Gravity with Perturba- 
tive Constraints (see also Ref. [24[ for an initial discussion 
of cosmology with perturbative constraints) . In contrast 
to the fundamental approach, we do not propose a basic 
theory but rather remain largely neutral on the details 
of how a particular set of observations are generated. In 
contrast to the phenomenological approach, our models 
have a wide range of applicability and one can, for ex- 
ample, ask whether the cosmological predictions that do 
not add extra degrees of freedom are consistent with ob- 
servations in the solar system or of compact objects. 

Our aim is to cover as large a class as possible of mod- 
ifications to gravity that do not require new degrees of 
freedom. Within our framework, it is possible in princi- 
ple to incorporate all such modifications that are covari- 
ant, metric (in the sense of obeying the Einstein equiv- 
alence principle), and conserve energy-momentum. The 
price we will pay for such generality is that we must re- 
main perturbatively close to General Relativity. While 
we discuss in detail the exact conditions for perturbative 
validity below, we caution here that this valid- ity does 
not require weak fields (as defined with reference to ei- 
ther the Minkowski or homogenous FRW metrics as the 



background) . 

Our method of generating such models is known in 
the literature as that of perturbative constraints [25[ ; it 
is an alternative way of associating, in a physically un- 
ambiguous fashion, a general action of the gravitational 
field with a field equation that is of second order. Devel- 
oped to deal with approximate local actions derived from 
fundamentally non-local theories, it has a wide range of 
applicability and does not rely on the existence of an 
underlying non-locality. 

The method of perturbative constraints we use here 
allows us to retain many of the important features of 
the action of the gravitational field. In particular, it al- 
lows us to retain the translational symmetries that, un- 
der Noether's principle, lead to energy-momentum con- 
servation, as well as enforcing diffeomorphism invariance 
- the invariance of physical law under changes of coor- 
dinates. It does so without introducing new degrees of 
freedom, while at the same time producing behavior that 
differs from the General Relativistic case, and allowing 
for a consistent set of predictions across a wide variety of 
spacetimes. 



II. GRAVITY WITH PERTURBATIVE 
CONSTRAINTS 

In this paper, we specialize to the case of a flat, ho- 
mogenous FRW metric, 



ds 



-dt 2 + a 2 (t) [dr 2 +r 2 (d6 2 + 1 



(1) 



It is important to emphasize here, however, that because 
we will be modifying gravity in ways that affect cosmo- 
logical expansion, standard relationships between, e.g., 
the energy density of various components and the rate 
of expansion that hold in a flat Universe within General 
Relativity will not necessarily hold here. 

We start with an action for the gravitational field that 
is a general function f(R) of the Ricci scalar curvature. 
The field equation for this general theory is 

f'R^ - ~W + g^n f - v„v„/' - torGT^ , (2) 

where f = §i- For a na ^ FRW spacetime, the Ricci 
scalar curvature is 



R = 6 



(3) 



The above field equation obeys the Bianchi identity and 
hence V^T^ V = 0, which results in the usual conservation 
law 



Pi + 3~{pi + Pi) = 
a 



(4) 



Here pi and Pi are the energy density and pressure, re- 
spectively of the i th constituent of the Universe. In this 
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paper, we will consider Lagrangian densities of the form 



f(R) = R-2A 



2(n+l) 



R" 



(5) 



where n and p are free parameters and the last term in 
the action is considered only at the perturbative level. 
The case A = 0, n = 1 was first discussed in Ref. [24j. 

Henceforth we shall adopt geometerized units where 
8itG = 1. Following standard procedure, we evaluate the 
(/!, v) = (0, 0) component of the field equation as well as 
its trace and obtain 



,d 1 ,, • d 

3/'- + -/ + 3/"R~=p 
a 2 a 



(6) 



and 



Rf - 2/ - 3 



f"R + R 2 f" + 3-Rf 
a 



= (-p + 3P). (7) 



In this last equation, p and P denote the total density 
and pressure of all constituents of the Universe. 

When the function f(R) is non linear in the Ricci scalar 
curvature, Eq. [6] is not a first-order ordinary differential 
equation in time, as is the standard Fricdmann equation. 
Indeed, as discussed earlier, a non-linear Lagrangian in- 
troduces additional degrees of freedom to the gravita- 
tional field. The key feature of gravity with perturbative 
constraints is a consistent way of treating these addi- 
tional degrees of freedom as fixed and not allowing them 
to evolve independently [H HE Hi] . 

In practice, we combine Eq. [6] with the conservation 
law of Eq. 21 and obtain a single differential equation for 
the scale factor that we formally denote by 



G(a, a) + p T{d, a, a, a) = 
Solving this equation with p 4 = 0, i.e., 
S(d(°>,a (0) )=0. 



(8) 



(9) 



and requiring the scale factor to be unity at the present 
epoch, we obtain the General Relativistic solution, which 
we denote by cS°\i). We then use this solution to evalu- 
ate the derivatives that correspond to all the new degrees 
of freedom in Eq. {81), i.e., 



e(d,a)+/i 4 ^(a (0) ,« (0) ,« (0) ,« (0) )=0 



(10) 



This last equation is again an ordinary differential equa- 
tion that is of first order in time and we solve it requiring 
the scale factor to be unity in the present epoch. 

Inspection of the field equations or of the Lagrangian 
action shows that we can expand the final solution for 
the scale factor in orders of the parameter ^ 2 ( n+1 ) 
where Rq — p m ,o + 4A is the unperturbed Ricci scalar 
curvature of the Universe in the present epoch. For the 
case with n = 1, we can write the solution for the scale 
factor in the form 



o(t) 



= „(°) 



(*) 



1 



(p m .o + 4A) 2 



,(4) 



(0 



(p 



m,0 



4A) 4 



,( 8 ) 



(*) 



.(11) 



Note that, although the modified Friedmann equation 
includes only terms of order /i 4 , its solution has higher- 
order terms because of the non-linearity of the differential 
equation. Consistency requires that terms of order p, 8 
and higher are negligible and can, therefore, be dropped. 

When discussing our results, we will occasionally ex- 
press them in terms of the contributions of matter and 
of the cosmological constant to the energy density of the 
Universe, which are given by 



and 



3H 2 
A 



3H 2 ' 



(12) 



(13) 



respectively, where H = (a /a). Finally, in order to com- 
pare our results to the usual parametrization of the dark 
energy equation of state, we will also write the Friedmann 
equation in the XCDM form [27j as 



(14) 



3 \-J = Pm + PX , 

and the equation for acceleration as 



6- = -p m - (1 + 3w x )px , 
a 



(15) 



where px is the energy density of an equivalent "dark en- 
ergy" component with an equation of state Px = wxPx- 
Note that this term includes contributions from both the 
cosmological constant and by the perturbative term that 
modifies the Einstein equations. 



A. Case 1: A = 0, n = 1 

This situation corresponds to a Universe with zero cos- 
mological constant and f(R) = R + p 4 /R. Since we are 
interested in deviations in the matter dominated regime, 
we will only consider matter with P = 0. In this case, 
the (/x, v) = (0, 0) component of the field equation and 
its trace become 



3/x 4 



R 

g" 



(16) 



and 



R 



3/i 4 
R 2 




-Si 



Pm 



(17) 
respectively. 

When p = 0, the equation for the scale factor describes 
the familiar matter-dominated Einstein cosmology, for 
which 



,(0) 



(*) 



y^Spmfit 



(18) 
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As a result 



Pm,0 

3 (aM) 



(19) 



and 



R<® = 6 



j(0) 



Pm,0 
ffl(0)V 



(20) 



Now to order [i 4 Eq. (fl~6| gives 



(i?(°)) 2 I 6 o(°) flW 

(4) Pm,Q 

fa(0)V 



(21) 



Solving for (t) using Eq. (JT9J) and ([20} the contribu- 
tion to order \i 4 is 



40 



(22) 



We then perform the same expansion up to order /i 8 . The 
contribution is 



»«(t) = 



18819p 4 ra ^ 
1600 



(23) 



As a boundary condition, we require that at the present 
epoch, which we denote by to, the scale factor is unity. 
Thus, as we evaluate a(t) to higher order in fi, we get 
corrections to the value of to from General Relativity, 



i.e. 



to 



l _ 3 n 



& Pm,0 



115 n 8 

2 Pm,0 



(24) 



As discussed earlier, it is practically impossible to 
check the convergence of our method without knowing 
the form of the Lagrangian action to higher order in the 
parameter /i. We can provide, however, a necessary con- 
dition for convergence in calculating an observable quan- 
tity, by requiring that the correction to the value of this 
observable calculated to order /x 8 is not larger than the 
correction calculated to order /x 4 . 

As an example, we calculate the value of the Hubble 
constant at the present epoch, H = H{t — t ), as 



Requiring the correction to order /i s in the expression 
for the Hubble constant to be smaller than the correc- 
tion to order /j, 4 results in the bound p 4 /p m o < 6/1017, 
whereas the same requirement for the deceleration pa- 
rameter leads to /x 4 /p m < 36/12312. Both of these con- 
straints at the present epoch are very similar and require 
/i 4 to have very small values, which cannot lead to an 
accelerated expansion within the limits of perturbative 
validity. 

An additional problem of using a cosmology with n = 1 
and A = to account for the observations can be seen by 
writing the Friedmann equation to order [i 4 (as in [24[), 
i.e., 



6p 4 
Pm 



(27) 



Transforming this equation in the XCDM form requires 
that 



Px 



Pm 



(28) 



Clearly the matter-dominated perturbative gravity the- 
ory to order /i behaves as a dark energy fluid with tux = 
—2 which is inconsistent with the combined WMAP and 
BAO data [|. 



B. Case 2: A ^ 0, n = 1 

We have found that, while the case with n — 1 and 
A = may generate accelerating solutions, these solu- 
tions lie outside the range where the perturbative expan- 
sion to the action is valid. Intuitively, this comes from 
the fact that by the time acceleration has begun, the mat- 
ter density is low enough that the zeroth-order relation, 
Ro = Pm.o, implies that the contribution of the addi- 
tional term, p 4 /R, is too large to be treated as a small 
perturbation to the standard gravitational action. This 
is not the case, however, in a Universe with a non-zero 
cosmological constant, as we will discuss below. 

In the case of a Universe with a non-zero cosmologi- 
cal constant, the (/i, v) — (0,0) component of the field 
equation and the trace are 



3/x 4 



R 
6" 



Pm + A , (29) 



H 



Pmfl 



Pm,0 



1017 / 

2 Pm.O 



(25) 



Similarly we evaluate the deceleration parameter go 
—aa/a 2 at the present epoch as 



'7n 



1 - 36 



P 



n 2 

Pm,0 



12312 



n 4 

Pm,0 



(26) 



and 



R 



3/x 4 




"Da 



Pm + 4A , 

(30) 
respectively. 

We will again seek a perturbative solution to this equa- 
tion of the form (jlip . When /i — 0, the Friedmann equa- 
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tion is the same as for ACDM, i.e 
a<°>(t) = 



(31) 



The time at the present epoch, i.e., the one for which the 
ACDM solution gives a scale factor of unity, is 



to 



sinh 




(32) 



At orders /i 4 and fi & , the solutions of the equation for 
the scale factor, a^(t) and a^(t), respectively, as well 
as the time to that corresponds to the present epoch can 
be found in closed form but are too long to be displayed 
here (we present the solution to order /i 4 in Appendix A) . 
Instead, we write the Fricdmann equation to order \i 



3A-6(/> m +A) 



4A 



(33) 



and the acceleration equation to the same order 

,4 



6- 

a 



-Pro + 2A 



fl- 



ip™ + 4A)2 
6 A + 24 ( Pm + A) 

-54(p m + A) 



Pn 



Pm 



Pm + 4A 

2 



4A 



(34) 



These two equations allow us to calculate to order p 4 two 
observable quantities at the present epoch, as a function 
of the parameters /i 4 and A: the Hubble constant, or 
equivalently the parameter f2 mj o, and the deceleration 
parameter go. We show the result of this calculation in 
Fig.Q]and[2]respectively, together with the regions where 
the calculation of each observable is not perturbatively 
valid. 

In both figures, the horizontal axis, i.e., /i 4 /(p m ,o + 
4A 2 ) = 0, corresponds to a General Relativistic Universe 
with a cosmological constant. The WMAP results, inter- 
preted within General Relativity, correspond to a Uni- 
verse with A/ p m , — 2.65 Keeping this value constant 
and moving towards positive values of /u 4 /(p m ,o + 4A 2 ), 
the predicted value of f2 mi o increases, because the per- 
turbative term in the modified Friedman equation con- 
tributes negatively to the rate of expansion of the uni- 
verse, and the deceleration parameter decreases for the 
same reason. The opposite is true when the parameter 
p 4 /(Pm,o + 4A 2 ) becomes increasingly negative. 

The hatched areas in both figures show the regions of 
the parameter space where the corrections to each observ- 
able calculated to order /i 8 are not negligible and hence 
the solution is no longer of the same order as the field 
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FIG. 1: Contours of constant values of the parameter f2 m 
at the present epoch, for different values of the cosmological 
constant A and the parameter /i 4 of the action of the gravi- 
tational field. The normalization constant p m ,o is the density 
of matter at the present epoch. The hatched area shows the 
region of the parameter space where the calculation of f2 m is 
not perturbatively valid. 
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FIG. 2: Same as Figure [T] for the deceleration parameter qo 
at the present epoch. Only negative values of qo are shown, 
which correspond to Universes that are currently accelerating. 



equation. This is the criterion we have discussed earlier 
in determining the perturbative validity of our calcula- 
tions. The particular shapes of the hatched regions are 
determined by the dependence of the perturbative terms 
in Eqs. [33] and [3U both of which can be negative, posi- 
tive, or even zero for different values of the ratio A/p mj0 - 
In both figures, however, the excluded regions lie far from 
the parameters of the General Relativistic Universe that 
are consistent with the WMAP data, suggesting that 
small potential deviations from the General Relativistic 
predictions can be modeled successfuly within our frame- 
work. 

Finally, we investigate the parameter u>x at the present 
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FIG. 3: Same as FigureQ]for the equivalent equation-of-state 
parameter wx of a XCDM model with the same cosmological 
expansion at the present epoch. 



As with Case 1, we can write the Friedmann equation 
to order /i 2 (™ +1 ) and obtain 



31* 

a 



2(n+l) 



2pS 



[6n(n + 1) + n- 1] 



(37) 



At the same order the acceleration equation is given by 
(2 + 3n)^ 2(,l+1) 



6- = -p„ 
a 



2p n 

r m 



[6n(n+l) + n-l] . (38) 



If we write the Friedmann equations in the flat XCDM 
form, Eq. (fFi]) , in the present epoch, we obtain for the 
equivalent equation-of-state parameter, 



w x 



-(n + 1) 



(39) 



As a result, increasing the value of the exponent n makes 
the equation-of-state parameter wx more negative. 



epoch for an equivalent XCDM cosmology that would re- 
sult in the same cosmological expansion rate and accel- 
eration as that predicted by the last two equations. We 
calculate the equivalent equation-of-state parameter by 
comparing the last two equations to those of the XCDM 
framework. The result is, to order /x 4 , 



w x 



-1 + 



12 



pm + A 



(p m + 4A) 2 V A 



- 18 



Pm + 4A 



Pn 



Pm + 4A 



(35) 



In Fig. [3l we show contours of constant values of the pa- 
rameter wx- The current observed limits —0.14 < wx < 
0.12 on the equation-of-state parameter H suggest that 
deviations from the ACDM model that do not introduce 
new degrees of freedom fall within the perturbative limits 
of our approach and cannot be ruled out a priori. 



C. Case 3: A = 0, n 1 

Finally, we examine a more general class of theories 
with n^l that are nevertheless matter dominated. We 
focus here on the case A = in order to isolate the ex- 
pected dependence of our previous results on the value 
of the parameter n. 

In this case, the (p, v) — (0, 0) field equation becomes 



3// 



2(1, 



+ n — h n(n + 1) 



(36) 

and its solution can be expanded again in terms of the 
perturbative parameter p 2 / Rq. The constrain on this 
parameter for perturbative validity becomes increasingly 
less restrictive as the exponent n decreases. 



III. CONCLUSIONS 

We have presented a formalism - Gravity with Pertur- 
bative Constraints - that describes deviations from Gen- 
eral Relativistic predictions at cosmological scales with- 
out introducing new degrees of freedom. Our formalism 
preserves diffeomorphism invariance, the Einstein equiva- 
lence principle, and energy-momentum conservation. Be- 
cause it is based on modifications of the Lagrangian of 
the gravitational field, it allows us to link a wide variety 
of observations - from compact objects and solar system 
phenomena, to non-linear clustering and linear structure 
formation - in order to test for the consistency of assum- 
ing no new gravitational degrees of freedom. 

Our method is quite different from the Palatini for- 
mulation of 1/R gravity discussed by, for example, 
Ref. 1331. and expanded on as "modified source gravity" 
in Ref. . Since both our methods and those related to 
Palatini formulations, however, eliminate the extra de- 
gree of freedom associated with modifying General Rel- 
ativity it is worth comparing some of their features. In 
particular, because we are able to treat our system per- 
turbatively - and, explicitly, because our predictions for 
densities p > p CI1 t do not demand a particular behavior 
when p <C p cr it - we do not run into the issues found 
by Ref. [l4| , where the microscopic discreteness expected 
for ordinary baryons and dark matter particles no longers 
"averages out" correctly to return the homogenous Fried- 
mann equations. 

Put explicitly, Ref. [l4[ noted first that the microscopic 
discreteness of matter did not strongly affect the averaged 
metric in General Relativity: i.e., when taking the mat- 
ter field to be composed of discrete particles (small but 
still larger than their Schwartzschild radius) the gravita- 
tional fields relative to the background are still small and 
g^ LV can be consistently be replaced by {g^ v ) plus a small, 
linear perturbation. This behavior was then contrasted 
with the Palatini 1/R formulation, where the $ field, a 
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function of the matter density that appears as a modifica- 
tion to the Einstein equation, has strong variations from 
point to point since the behavior of <E> is strongly non- 
linear in the range p £ [0, p C rit] regardless of the length 
scale of the perturbation. 

In the case of gravity with perturbative constraints, 
however, describing the behavior of the equations of mo- 
tion when p <C jOcrit may involve increasingly higher or- 
der terms in the approximate equations of motion, and 
- likely for many, though perhaps not all, models where 
the first terms in the equations of motion are generated 
by l/R n cases discussed in this paper - eventually goes 
outside the radius of convergence of the series. That our 
theory can not in all cases make predictions for atomic- 
scale physics is a feature of the perturbative approach; in 
contrast to Palatini 1/R, the behavior at p > p cr j t decou- 
ples from these more difficult questions and no particular 
atomic-scale behavior is required by us for consistency 
with larger scales. 

In this first approach to the problem, we chose to look 
for modifications of the Einstein-Hilbert action that in- 
volve the addition of terms proportional to an inverse 
power of the Ricci scalar curvature. We were lead to 
this choice by our goal to describe potential deviations 
at cosmological scales, without affecting significantly the 
behavior of gravity in the solar system. One might ob- 
ject that our particular choice of parametrization - a two 
dimensional space, (A, n) - too-drastically narrows down 
the possible models considered. In a sense, this is a prob- 
lem with any attempt to cover a space of functions with 
a finite number of parameters. 

Our approach, however, follows the general direction 
of testing for modifications of gravity in the solar sys- 
tem and in other astrophysical systems by using the 
Parametrized Post-Newtonian (PPN) framework [28| • In 
a similar fashion to our work, the PPN framework pro- 
vided a way to link together different observations in var- 
ious astrophysical settings in a consistent and physically 
rigorous way. While the PPN model itself required a cer- 
tain degree of arbitrariness - the functional forms of the 
potentials - it turned out that the original choices of these 
functional forms [29| were sufficient to cover nearly every 
theory proposed in the following thirty-seven years. Only 
in a few cases have failures of PPN been due to an overly 
restrictive choice of parameters (see, e.g. 30]). The use- 
fulness of the framework we propose here will depend, of 
course, on its ability to capture the behavior of plausible 
gravitational theories that go beyond General Relativity. 

The particular set of terms that we included in the La- 
grangian of the gravitational field allow us also to study 
in a more systematic way other phenomenological cos- 
mologies that produce accelerating expansion. For ex- 
ample, the set of parameters discussed in Sec. Ill CI i.e., 
A = and n^l, provide a homogenous cosmology iden- 
tical to the Cardassian model of Ref. [HI]. Within our 
formalism, however, we can also make predictions for 
the formation of structure in this cosmology, by taking 



explicitly into account the modifications to the Poisson 
equation. 

Our formalism is also capable of describing the cos- 
mological expansion histories that are generated by the 
phenomenological XCDM model [3l|. This is demon- 
strated in Fig. [31 for a flat universe, where the values of 
the equivalent equation-of-state parameter w are shown 
as a function of the perturbative parameter p 4 , the mat- 
ter density, and the value of the cosmological constant. 
In a similar manner, our framework can also describe 
the cosmological expansion generated by the model sug- 
gested by the Dark Energy Task Force [10| , in which the 
equation-of-state parameter is allowed to change in time 
as wq+(1— a)w a , since the values of the equivalent param- 
eter uu calculated here depend indeed on cosmic time (cf. 
Eq. [35]) . In contrast to the phenomenological model, how- 
ever, our framework offers additional insight as to how 
the resulting gravity modifications agree with structure 
formation observations, or with "fifth-force" constraints 
from other arenas. 

We will study the behavior of the cosmological expan- 
sion predicted by gravity with perturbative constraints 
for a Universe with a finite spatial curvature as well 
as compare these predictions directly to observations in 
forthcoming papers. 
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IV. APPENDIX: PERTURBATIVE SOLUTION 
TO 4TH ORDER FOR A # 

As discussed in the main text, perturbations to ACDM 
at order p 4 has closed-form solution. In particular: 



2tan~ 1 


V^3 tanh ^ 




- 9V3t 


96^3 tanh 







17 - 23 cosh(V3Ai) + 8 cosh(2V3Ai) 



48 1 - 2cosh(V3At) 
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